Excited nucleons with chirally improved fermions 
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We study positive and negative parity nucleons on the lattice using the chirally improved lattice 
Dirac operator. Our analysis is based on a set of three operators Xi with the nucleon quantum 
numbers but in different representations of the chiral group and with different diquark content. We 
use a variational method to separate ground state and excited states and determine the mixing 
coefficients for the optimal nucleon operators in terms of the x*- We clearly identify the negative 
parity resonances A^(f535) and A'^(i650) and their masses agree well with experimental data. The 
mass of the observed excited positive parity state is too high to be interpreted as the Roper state. 
Our results for the mixing coefficients indicate that chiral symmetry is important for A'^(i535) and 
A'^(i650) states. We confront our data for the mixing coefficients with quark models and provide 
insights into the physics of the nucleon system and the nature of strong decays. 
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I. INTRODUCTION 

Lattice QCD has the potential to answer some of the 
most important questions about the proper effective de- 
grees of freedom in the nucleon and its excitations. In 
particular, one is able to arbitrarily change the mass 
of quarks and hence to study the evolution of nucleon 
properties such as its mass and excitation spectrum as 
a function of the current quark mass. This is interest- 
ing because we know the physical picture which drives 
hadrons made of heavy quarks. In this case the hadrons 
can be well approximated by the non-relativistic QCD 
picture, where the most important ingredients are the 
linear -I- Coulomb potential as well as a hyperfine correc- 
tion that originates from the perturbative color-magnetic 
interaction between spins of heavy quarks. 

The most prominent implication of this picture is the 
structure of the heavy quarkonium spectrum. The color- 
magnetic spin-spin force induces a splitting between the 
^Si and ^5*0 ground states. The linear + Coulomb inter- 
action drives the pattern of radially and orbitally excited 
states: The first orbital excitation which carries L = 1 
orbital angular momentum and has a parity which is op- 
posite to the given ground state, lies always below the 
first radial excitation, which has the same quantum num- 
bers as the ground state. Such a pattern is quantitatively 
understood within non-relativistic lattice QCD as well as 
the potential approach (for a review and references see 
e.g. Ref. X\). 

On the other hand the ordering of the lowest positive 
and negative parity excitations of the nucleon is just op- 
posite. The first radial excitation of the nucleon, the 
Roper resonance 1/2+, Ar(1440), lies about 100 MeV be- 
low the first orbital excitation, l/2~, iV(1535). A similar 
pattern is observed in the A spectrum, while the opposite 
ordering of the lowest positive-negative parity J = 1/2 



states takes place in the A spectrum. It appears to be im- 
possible to explain this within a linear confinement plus 
perturbative gluon exchange correction picture. 

The whole experience of nuclear physics as well as nu- 
merous data on the chiral structure of the nucleon teach 
us that the consequences of spontaneous chiral symme- 
try breaking should be very important for the nucleon. 
A possible - certainly oversimplified - picture is that in 
the low energy regime the nucleon can be approximated 
as a system of three quasi-particles with dynamical mass 
induced by spontaneous chiral symmetry breaking (con- 
stituent quarks) which are confined. Such a picture is the 
most straightforward realization of the SU{(!>) symmetry. 
The latter, in turn, can be proved to be an exact symme- 
try of QCD for ground state baryons in the large Nc limit 
yf. Then the question arises, what are the other residual 
interactions between the constituent quarks which break 
SU{6) symmetry? The empirical ordering of positive and 
negative parity states is opposite in mesons and baryons 
which indicates that some physical aspects are different 
in these systems. An anomalous ordering in baryons was 
explained to be possibly due to a significant contribution 
of the residual interaction mediated by the chiral (Gold- 
stone boson) field Q . A meson-like exchange interaction 
between valence quarks is possible within the quenched 
approximation in baryons, and not possible in mesons. 

Attempts to address these issues within the lattice for- 
mulation have begun only recently 0, IE IE S IE IS ^| ■ 
Mainly the two interpolating fields xi ^^nd X2 (for their 
definition see below) have been used in order to extract 
the ground state nucleon as well as the lowest excited 
states of both parities. The operator xi is known to cou- 
ple well to the nucleon ground state and has been used in 
lattice calculations for more than 20 years. This opera- 
tor also couples to the = l/2~ (negative parity) state. 
And indeed, the lowest 1/2^ state was established and 
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has a tendency to approach the physical iV"(1535) state 
for smaU quark mass. The authors of see no sphtting 
between the lowest and the next excited negative par- 
ity state = l/2~ Af(1650), whereas such splitting is 
noticed in 0,0- 

The first excited positive parity iV(1440} Roper state 
has not been identified in lattice studies 0,|a,0>l3 • Initial 
hopes that it would have a reasonable overlap with X2 
have been abandoned since the plateaus seen from the 
X2 correlator come out consistently too high (close to 2 
GeV) to be the Roper state. 

Recently |^ the Roper state has been observed and 
reordering of the levels, such that the excited positive 
parity state falls below the negative parity states taking 
place at pseudoscalar masses somewhere between 200 - 
500 MeV, was obtained in '3, % 1^ using the xi inter- 
polator. Both articles rely crucially on advanced fitting 
techniques, i.e. Bayesian priors. 

Here we tackle the analysis of the nucleon system com- 
bining different powerful strategies. In addition to the 
operators Xi,X2 we use a third operator Xa and com- 
pute all cross correlations. This correlation matrix is 
then analyzed using the variational method [ILj which 
allows for a separation of ground and excited states. Us- 
ing more operators in the variational method allows the 
system to better select those combinations of interpola- 
tors which have optimal overlap with the physical states. 
The possibility of a reordering of the levels at small quark 
masses suggests to use a Dirac operator with good chiral- 
ity which allows one to go to small quark masses. Here we 
work with the chirally improved Dirac operator [l^ . We 
perform a careful analysis of finite volume effects which is 
crucial in a study of excited nucleons with their extended 
wavefunctions. 

From the variational technique one obtains also im- 
portant information on the optimal operators which have 
maximal overlap with the physical states. This is inter- 
esting not only from a technical point of view but it also 
provides important insight into the physics of the nu- 
cleon system. We compute the mixing coefficients which 
relate our operators Xi to the optimal operators, i.e. the 
coefficients which determine the operator content of the 
physical states. These mixing coefficients are then stud- 
ied as a function of the quark mass and confronted with 
information from quark models. We discuss the physical 
implications of this comparison. 

Our article is organized as follows: In the next sec- 
tion we collect information on the lattice calculation such 
as simulation details, the operators used, the variational 
technique and our fitting methods. In Section 3 we dis- 
cuss the diquark content of our nucleon interpolators and 
their chiral properties. This is followed by two sections 
where we discuss our raw data and then present the re- 
sults for the mass spectrum. In Section [VII we determine 
the mixing coefficients, discuss their physical implications 
in Section IVIII and we end with conclusions in Section 
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II. LATTICE TECHNOLOGY 

For our quenched calculation we use uncorrelated con- 
figurations from the Liischer-Weisz gauge action ^3 ■ We 
work on 16^ x 32 and 12^ x 24 lattices using 100 config- 
urations for both volumes and gauge coupling constant 
Plw = 7.9. The lattice spacing is a = 0.15 fm as de- 
termined from the Sommer parameter in We re- 
mark that the scale from the Sommer parameter agrees 
well with a determination from hadronic input |l5j | . This 
gives rise to spatial extensions of 2.4 fm, respectively 1.8 
fm. Throughout this paper we use this scale where we 
give mass values in physical units. 

The gauge field configurations were treated with one 
step of hypercubic blocking 16]. We use the recently 
developed Chirally Improved (CI) Dirac operator [l^ . 
The CI operator is a systematic expansion of a solu- 
tion of the Ginsparg- Wilson relation and has very 
good chiral properties. In |15i] the residual quark mass 
for the Plw = 7.9 ensembles used here was computed 
from the chiral extrapolation of the pion mass and found 
to be nires = -0.0020(5)a"i = -2.7(7) MeV. It was 
demonstrated 0,0] that the CI operator allows simula- 
tions with pseudoscalar-mass to vector-mass ratios down 
to mps/my = 0.28 at relatively small cost and that it 
has good scaling properties. It was found 0| that for 
the nucleon the scaling dependence on a is of the size 
of the statistical error. For the states we focus on here, 
i.e. the excited positive parity state and the negative 
parity states, we have larger errors and the effects of fi- 
nite a are negligible. We use Jacobi smeared sources 
1^ with a width of about 0.7 fm. For the nucleon we 
find that the signal at small quark masses is enhanced 
when smearing also the sink. This allows to go deeper 
to the chiral limit, while at the larger quark masses we 
find exact agreement between the data from smeared and 
point-like sinks. Thus for the nucleon we use smeared 
sinks for all quark masses. All correlators were projected 
to zero momentum. The inversion of the Dirac opera- 
tor was done with the BiCGstab multi-mass solver j23| 
and we use several different values of the quark mass in 
the range between aniq = 0.013 (mpg — 270 MeV) and 
airiq = 0.2 {mps = 866 MeV). A more detailed account 
of our setting is given in 15]. 

For the extraction of nucleon data we implement the 
following set of operators: 

Xi{x) = eabc[ua{x)Cj5di,{x)]uc{x), (1) 
X2{x) = eabc[ul{x)Cdb{x)] j5Uc{x) , (2) 
X3(a;) = ieabc[uaix)C jf^^5db{x)]uc{x) . (3) 

C is the charge conjugation matrix and T denotes trans- 
position. The color indices are a, b, c, while for the Dirac 
indices we use matrix/vector notation. For the interpo- 
lator X3 we use only the time-like component. All three 
operators have the quantum numbers of the nucleon but 
correspond to different representations of the chiral group 
and contain different diquark fields in the brackets in Eqs. 
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(0 - (O. Therefore one would expect that they couple 
differently to the given physical state. In the next sec- 
tion we will discuss in detail these representations, the 
quantum numbers of the diquark fields and their physi- 
cal implications. 

We use the variational technique compute all 

cross correlations of our operators Xi listed in Eq. ^ - 
®: 



(4) 



We explicitly project our correlators to definite parity 
by inserting the projectors = ^[1 ± 74] and the 
trace in Eq. |0J is taken in Dirac space. Note that 4 
is the euclidean time direction and 74 the correspond- 
ing element of the Clifford algebra. In the following 
we often omit the superscript ± for brevity. In or- 
der to extract the ground state and excited levels we 
first normalize the correlation matrix M and construct 
D{t) = M-i/2(io)M(t)M-i/2(to), where the position to 
of the time slice chosen for the normalization can be used 
to optimize the signal. For all results presented here we 
have used to = 1, which in our conventions corresponds to 
propagation distance 0. We discuss this choice in Sectlvl 
Note that M is a hermitian, positive definite matrix 
such that its inverse square root exists. We remark that 
this normalization is designed to reduce the error but a 
separation of ground and excited states can be obtained 
also without this normalization. Actually we find that 
in our case we could work almost equally well with the 
unnormalized matrix M. The second step is a diagonal- 
ization of the hermitian matrix D{t) giving rise to three 
real eigenvalues Xi(t),i — 1,2,3. Applying the lemma 
proven in 11] one finds 



A,(t) = aexp[-{t~to)E,] 

xh+ofexp [-{t~to)AE^] 



(5) 



where AEi is the difference of Ei to the first omitted 
energy level, i.e. AEi = \Ei — £'4]. We remark that the 
form quoted in Eq. Q is the expression for the positive 
parity states when using the matrix M+, projected to 
positive parity. For fitting the negative parity states one 
has to time-reverse the correlators, i.e. replace thyT — t, 
where T denotes the temporal extension of the lattice. 
The eigenvalues from M~ are related to the eigenvalues 
from Af+ by a minus sign and time reversal. We av- 
erage our results from M+ and M~ after the necessary 
transformation to increase the statistics. 

Eq. (O) makes clear the strength of the variational 
method: Each eigenvalue corresponds to a different en- 
ergy level which dominates its exponential decay. For 
a large enough set of operators this allows for a clean 
separation of ground state and excited levels and we can 
apply a simple stable 2-parameter (Ci, Ei) fit to deter- 
mine the corresponding energy. In principle an N x N 
correlation matrix is sufficient for determining the ground 



state and the first A'^ — 1 excited states as long as the op- 
erators are linearly independent and each of the physical 
states has overlap with at least one operator. However, 
using more operators improves the results. The "perfect 
operator" has infinitely many different lattice operators 
(point split, different levels of smearing etc.). Some have 
large coefficients and others don't. Selecting those with 
large overlap is a highly non-trivial task and requires in 
principle a knowledge of the wave function. 

For a somewhat simpler problem - the application of 
the correlation matrix technique to a scattering problem 
- where there is a natural ordering of the operators (rel- 
ative momenta of the scattering particles) a criterion to 
determine how many operators are needed was given in 

m 

Our functional uses the full correlation matrix. The 
fit ranges were determined from effective mass plateaus 
and we find values of x^/d.o./. ranging from 0.5 to 1.5. 
Slight increases or decreases of the fit range leave our fit 
results and the /d.o.f. essentially invariant. We found 
that in the negative parity sector and for the excited 
positive parity state the quality of the data decreases for 
small quark masses and we could not maintain the stan- 
dards of our fitting procedure. We do not give our results 
for these cases. The errors we quote were computed using 
jackknife subensembles. 



III. DIQUARK CONTENT OF NUCLEON 
INTERPOLATORS AND THEIR CHIRAL 
PROPERTIES 

Let us now discuss in more detail the structure of our 
nucleon operators Xi defined in Eqs. CD " ®- The fields 
Xi and X2 create from the vacuum both, positive and 
negative parity states, with I, J = 1/2, 1/2. Linear com- 
binations of these fields have intensely been used in QCD 
sum rule calculations . The field xi contains a scalar 
diquark [s^ I, = 0, O"*" component, while the field X2 
has a pseudoscalar diquark component /, =0, 0^. 

It is instructive to review the chiral SU{2)l x SU{2)ji 
properties of these fields. Here and in the follow- 
ing (/i, /fl) denotes the irreducible representations of 
SU{2)l X SU{2)r with II and Ir being isospins of the 
left-handed and right-handed components. The quark 
field q{x) (with q = u or d) together with its parity part- 
ner 75(7(3;) forms a (0, 1/2) © (1/2,0) reducible represen- 
tation of the chiral group. This reducible representation 
is an irreducible representation of the full parity-chiral 
group SU{2)l X SU{2)ii x d, where the group Qcon- 
sists of the elements identity and space inversion |23l | . 

Consider first the scalar diquark field in This di- 
quark field is invariant with respect to both isospin as 
well as axial transformations in the isospin space and 
thus transforms as a scalar representation (0, 0) of the 
chiral group. The chiral transformation properties of 
the interpolator xi are then determined exclusively by 
the transformation properties of the third quark. Hence 
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the interpolator xi together with its parity partner 75 xi 
transform as (0, 1/2) © (1/2, 0), 

{Xi(a^), 75 Xi(x)} ~ (0, 1/2) © (1/2, 0) . (6) 

Similarly, the pseudoscalar diquark in X2 also is invariant 
with respect to chiral transformations. Again the trans- 
formation properties of X2 are solely determined by the 
transformation properties of the third quark field, 

{X2(a;), 75X2(x)} ~ (0, 1/2) © (1/2, 0) . (7) 

We stress that the fields {xi{x), 75X1(2^)} and {x2{x), 
75 X2ix)} belong to distinct irreducible representations of 
the parity-chiral group and they do not transform into 
each other under chiral transformations. Linear combi- 
nations of the fields xi and X2, as used in the QCD sum 
rule approach, form only reducible representations of the 
parity-chiral group ^3 ■ 

The transformation properties with respect to a flavor- 
singlet U{1)a transformation are quite different, because 
the flavor-singlet axial transformation mixes scalar and 
pseudoscalar diquarks. Only the specific linear combina- 
tions 

{(Xi(a;)+X2(a;)) , «75 (Xi(a^)+X2(a;))} , (8) 
{ (Xi(a;) - X2(a;)) , «75 (Xi(a;) - X2(a;)) } , (9) 

form distinct irreducible representations of U{1)a x Ci. 

There is no reason to expect that the interpolators xi 
and X2 are the only ones that couple to the nuclcon and 
its excited states. Hence a-priori one should expect that 
the nucleon and its excited states are actually some mix- 
tures of different irreducible representations. In order to 
check this aspect we also included our interpolator ^3. 

In X3 the diquark field in the brackets has isoscalar- 
vector quantum numbers I,J^ = 0, 1^ for the spatial 
components (/i = 1,2,3) of the 4-vector and isoscalar- 
scalar quantum numbers I, = 0, 0+ for the time com- 
ponent {ji = 4). Hence for /i = 4 this operator cou- 
ples only to J = 1/2,/ — 1/2 states, similar to Xi and 
X2, whereas for ^ = 1,2,3 it couples to both J = 1/2 
and J = 3/2 (of both isospins 7 = 1/2 and / = 3/2) 
baryons. Contrary to the isoscalar-scalar diquark in xi 
the diquark field in X3 transforms as (1/2, 1/2) under 
SU{2)l X SU{2)i^ and these transformation properties 
distinguish it from the diquark in xi- This also im- 
plies that X3 together with some other fields generates 
the (1, 1/2) © (1/2, 1) and (0, 1/2) © (1/2, 0) represen- 
tations of the chiral group, the latter is however distinct 
from the representations generated by xi and X2- 

IV. DISCUSSION OF THE RAW DATA 

Before presenting results for the masses let us briefly 
discuss our raw data. This discussion illuminates the 
roles of the different interpolators Xi and underlines the 



strength of the diagonalization method. We base this 
analysis on effective masses computed from either the in- 
dependent diagonal correlation functions (xiXi) or from 
the eigenvalues A^. 

In Fig. n we show in the left-hand-side (Ihs.) column 
plots of effective masses as obtained from the indepen- 
dent diagonal correlators (xiXi), (X2X2) and (xaXs) (top 
to bottom). The right-hand-side (rhs.) column shows ef- 
fective masses for the three eigenvalues Ai, A2 and A3 of 
the correlation matrix. The horizontal bars shown in the 
plots for Ai and A2 indicate our fit results and the fitting 
range. The characters in the plots label structures which 
we discuss below. All data are from the 16^ x 32 lattice 
at a rriq = 0.10. 

Let us first discuss the independent diagonal correla- 
tors (xiXi) shown in the Ihs. column. The positive parity 
states are extracted from the left part {t < 18) of the 
plots, while the negative parity states that couple to the 
same interpolator dominate the right part {t > 18) of 
the plots. One clearly sees a good plateau for the inter- 
polators xi and X3 iii the interval (5, 15). This plateau 
corresponds to the ground state nucleon (positions a and 
e). At small Euclidean times one can see contributions 
from the excited states of positive parity. If an interpola- 
tor couples strongly to the nucleon ground state (as it is 
the case for xi and xa) then the signal from the excited 
states is only weak and the extraction of the mass for the 
excited states becomes a complicated task. 

The signal obtained with the interpolator X2 is intrigu- 
ing. There is no plateau corresponding to the nucleon 
ground state. This fact has a natural explanation within 
a quark model and will be discussed below. However, 
for relatively large quark masses aruq > 0.08 one sees 
a plateau at smaller distances (3,7), i.e. position c. We 
remark that for the arriq — 0.10 shown in the figure the 
quality of this plateau is already very poor, but it is more 
pronounced at our larger quark masses. The mass corre- 
sponding to this plateau is higher than the mass of the 
negative parity state seen from the same interpolator. 
Though we do not know the fate of this excited positive 
parity state at small quark masses, it is unlikely to be 
the Roper (cf. our discussion in Section rvil)l . 

Let us now look at the signals for negative parity 
states from the diagonal correlators (xiXi)- three 
of them show plateaus for time intervals (25,30), i.e. 
positions 6, d, /. Similar to the plateaus for the ex- 
cited positive parity state these plateaus are more pro- 
nounced for larger quark masses and vanish entirely be- 
low a rriq — 0.06. It is interesting to note that the plateau 
in (X2X2) (position d) is somewhat higher than the other 
two plateaus (positions fc, /). This is consistent with an 
observation we present below where we will show that at 
relatively large quark masses the heavier negative par- 
ity state Af(1650) has a larger X2 component, while this 
component is smaller for the lighter 7V(1535). 

The different states which appear in the individual 
diagonal correlators are separated when analyzing the 
eigenvalues of the correlation matrix. The correspond- 
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diagonal correlators < X; X; > 



eigenvalues X. 
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FIG. 1: Effective mass plots for the 16"^ x 32 ensemble at amq — 0.10. The Ihs. column shows results from the independent 
diagonal correlators (xiXi)j (X2X2) and (xsXs) (loP to bottom) while the rhs. column is for the eigenvalues Xi of the correlation 
matrix. The horizontal bars show our fitting results and range. The characters in the plots label structures which we discuss 
in the text. The insert in the top right panel demonstrates the effect of turning off the smearing of the sink. 



ing effective mass plots are shown in the rhs. column 
of Fig. n We order the eigenvalues Xi{t) with respect 
to their size, such that for each value of t we have 
Xi{t) > X2{t) > X^it). The nucleon ground state gives 
rise to a well pronounced plateau in Ai in the time in- 
terval (5,15) (position g). For the first excited positive 
parity state there is a hint of a plateau in the range (4,8) 
(position i), but the quality is as poor as for the cor- 
relator (X2X2) alone. We remark that there is no sign 
of the excited states contribution seen in the Ai plot at 
small times. The quality of the plateaus for the two neg- 
ative parity states (positions h and j) is somewhat bet- 
ter. These plateaus become more pronounced for larger 
quark masses but vanish entirely for aniq < 0.06. As 
is obvious from the bottom right plot the third eigen- 
value A3 is highly contaminated with excited states, has 



large statistical errors, and no structure is visible in the 
corresponding effective mass plots. 

For our variational analysis the number of independent 
operators should be sufficiently large and was limited to 
three for practical reasons. However, we also check our 
results by analyzing also the (xi, X2)-subsector. It turns 
out that the resulting eigenvalues are close to those in 
Fig-d Actually xs seems to contribute to the same states 
as represented by the other two operators and no clear 
signal to further states is observed. This is further dis- 
cussed in Sect. I VIII but also seen in Fig. where A3 
shows no additional plateau signal. 

In the insert (window) in the upper right plot of Fig. 
^we show the effect of omitting the Jacobi smearing at 
the sink. All the correlators shown in Figs. Handle are 
obtained with both smeared source and sink. The Ja- 
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FIG. 2: The eigenvalues Xi{t) and X2{t) on a logarithmic 
scale. The data are for 16^ x 32 at arriq — 0.04). The full line 
in the Ai-plot (top) shows the fit for the nucleon mass. 



cobi smearing is designed to improve the overlap with the 
ground state and the insert demonstrates clearly, that the 
contribution of excited states increases when omitting the 
smearing of the sink. For the point-type sink one now ob- 
serves a stronger signal from excited states at small time. 
Note also that since this signal appears in the Ai plot, and 
the A2 plot still exhibits a plateau (which originates from 
the X2 correlator), we conclude that this signal comes 
from the xi interpolator and the corresponding excited 
state is not the same as seen with the X2 interpolator. 
We have tried to analyze the xi correlator with two- and 
multi-exponential fits in an attempt to extract this pos- 
itive parity excited state. However, the statistical errors 
are too large for a reliable extraction and we cannot draw 
firm conclusions about the ordering relative to the neg- 
ative parity state. Hence we cannot decide from such a 
fit about the energy of the lowest positive parity excited 
state seen with the smeared source and point sink. While 
we cannot conclude anything concrete whether this signal 
comes from the Roper resonance or not, the insert shows 
that optimizing the spatial structure of the nonlocal cor- 
relators (in order to possibly separate the signals from 
the nucleon and the Roper) and introducing more oper- 
ators is one of the possibilities for further lattice studies 
of the Roper state. 

Let us finally comment on an interesting structure 
which we observe both in the X2 correlator and in A2 
for small quark masses. In Fig. |21 we illustrate this ef- 
fect by comparing Xi{t) and X2{t) on a logarithmic scale 
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FIG. 3: Masses of positive and negative parity states in GeV 
and in lattice units as a function of the pseudoscalar mass 
squared in GeV^ (quark mass in lattice units on top scale). 
We show data from 16'^ x 32 lattices as filled symbols and 
compare them to data from 12^ x 24 lattices (crosses). The 
curves represent the chiral extrapolation of the nucleon and 
the two negative parity states. The open symbols show the 
experimental values. 



for the 16"^ x 32 ensemble at anig — 0.04. The top plot 
shows Xi{t) and the full line indicates where we fitted 
the nucleon mass. The behavior of the second eigenvalue 
X2{t) (bottom plot) shows a somewhat surprising behav- 
ior: It starts off with a large slope but near t = 7 turns 
into a nearly linear behavior with a much smaller slope, 
in particular its absolute value is smaller than the one 
for the nucleon. If one interprets this as the trace of a 
physical positive parity particle it would be lighter than 
the nucleon. As a consequence we cannot interpret this 
structure as a signal from a physical state, but likely it is 
related to some quenched artifact. We remark that also 
the time behavior of the corresponding eigenvector (see 
Sect. IVI|) shows a change near t = 7 indicating that the 
left half of the A2 (t) does not describe a proper physical 
state. 



V. RESULTS FOR THE MASS SPECTRUM 

Let us now come to the presentation of the nucleon 
masses which we extract from our correlators. In Fig. 
13 we show the baryon masses in lattice units as a func- 
tion of the quark mass. Only data points were included 
where the standards of our fits discussed in Section |n] 
could be maintained and where a credible signal in the 
effective mass plot is seen. They all come from the diag- 
onalization approach using <o = 1- We did try up to 3. 
The quality of the plateaus does not improve and the re- 
sults remain invariant. Note that we use Jacobi smeared 
sources and there we find that these effects are minimal. 
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TABLE I: Fit results of fits to Ai and A2. The nucleon ground state is obtained from smeared-smeared results, whereas the 
remaining states are from smeared-point calculations. The mass values are in lattice units, the fit interval is denoted by t. 
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This behavior changes for more extended sources. All fit 
results and parameters are listed in Table 

The lowest lying set of data corresponds to the ground 
state in the positive parity sector, i.e. the nucleon. This 
is followed by two negative parity states which upon 
simple chiral extrapolation we identify as A'^(1535) and 
A^(1650). For brevity we use this particle data book 
nomenclature in the following. Although the two neg- 
ative parity states are relatively close to each other we 
clearly separate them with the variational technique dis- 
cussed above. The fact that the two negative parity 
states become approximately degenerate towards larger 
quark masses is well understood in potential models for 
heavy quarks. Finally, the excited state with positive 
parity is found on top and we will refer to it as . While 
we can identify the nucleon down to rather small quark 



masses, there is no clear signal (effective mass plateau) 
for the negative parity states and for the N'^ at small 
quark masses. However, below we will show that all 
states, except for the excited positive parity state N^, 
extrapolate well to their physical masses. 

For a study of excited states an analysis of finite vol- 
ume effects is particularly important. Excited states have 
more extended wave functions than the ground state and 
thus suffer more easily from squeezing them into a too 
small volume. Such a squeezing typically increases the 
kinetic energy of the constituent quarks and so pushes 
up the mass of the state. In order to analyze the finite 
volume effects, in Fig.|31we compare our numbers for the 
baryon masses as a function of the bare quark mass rUq 
for 16'^ X 32 and 12^^ x 24 lattices and find little volume 
dependence. 
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Our nucleon data easily penetrate the low pion mass re- 
gion where level crossing of the negative and positive par- 
ity states has been claimed in [IjBiEl- However, we find 
that for the other states the quality of our data depletes 
with decreasing quark masses and we could not maintain 
the standards of our fits for arriq < 0.06. We analyzed 
our data also with other methods (multi-exponential fits, 
Bayesian priors, maximum entropy method) but did not 
obtain stable, convincing results for an excited positive 
parity state below quark masses of am, = 0.06 as 
well as for the other positive parity excited state ex- 
tracted form the xi correlator. Although |^ |^ use 
slightly larger lattices (3.0 fm as opposed to our 2.4 fm) 
we cannot blame finite size effects for the depletion of the 
quality of our data, since we obtain essentially unchanged 
results already on our 1.8 fm lattice. 

The masses of the three states nucleon, A^(1535) and 
7V(1650) are extrapolated to the chiral limit. For the nu- 
cleon this is done by a linear extrapolation of the square 
of the nucleon mass (oAfjv)^ as a function of the quark 
mass amq. For the negative parity states we extrapolate 
linearly the mass as a function of the quark mass. Since 
the residual quark mass is so small (a rrires = —0.0020(5)) 
it was neglected in our chiral extrapolation. We used 
fully correlated fits to take into account the fact that the 
data for different quark masses were calculated on the 
same ensemble of configurations and thus are not inde- 
pendent. All three fits have a /d.o.f. less than 1. The 
linear extrapolation neglects chiral non-analyticities. We 
remark that we also experimented with functional forms 
suggested by quenched chiral perturbation theory |24j . 
However, our numerical data are not sufficiently precise 
for a clear distinction between the simple fits described 
above and the form from quenched chiral perturbation 
theory. Thus we decided to stay with the simpler 2- 
parameter fits. 

To set the scale we use either the p-mass (from [Tsj l') 
or give our results as ratios. In Table ^1 we list our fi- 
nal numbers for the masses in the chiral limit. We find 
that the nucleon mass matches its physical value reason- 
ably well and the two negative parity states are about 
8 % and 9 % too high. Given the fact that this is a 
quenched calculation the agreement of our results with 
experimental data is good. In order to estimate the sys- 
tematic error from the scale setting we compare the scale 
from the p-mass to the scale from the Sommer parameter 
(ap = 0.154(8) fm, ar„ = 0.148(2) fm). The second error 
given in the dimensionful results in the first line of Table 
nil was computed from the difference of the two scales. 

Like other authors 0, B Sll H II3- we also ob- 
serve a smooth increase of the N - A^(1535) splitting 
towards the chiral limit. At large current quark masses 
the hadrons may be described as a system with orbital 
motion in a color-electric confining field and the split- 
ting is due to this orbital excitation. The splitting slowly 
increases towards the chiral limit implying the increas- 
ing importance of another mechanism (spontaneous chi- 
ral symmetry breaking) in addition to confinement. This 



TABLE II; Comparison of our 16^ x 32 results in the chiral 
limit to experimental data. Masses are given in MeV, the 
ratios are dimensionless numbers. 
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data 
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1.63(2) 
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is consistent with the chiral constituent quark model [3j 
where an appreciable part of the splitting is related to 
the flavor-spin interaction between valence quarks. 

VI. THE OPERATOR CONTENT OF 
PHYSICAL STATES 

Not only the eigenvalues of the correlation matrix M 
in Eq. but also the corresponding eigenvectors pro- 
vide interesting information. In particular the entries of 
the eigenvectors give the mixing coefficients for the op- 
timal operators Xi which have maximal overlap with the 
physical states in terms of the original basis Xi- Let us 
be a little bit more explicit. The unnormalized matrix 
Mi,j {t) of Eq. Q is diagonalized by a unitary matrix U 
which is built from the eigenvectors 6*^'^ i.e. 

U M{t) C/t =diag(Ai(t),A2(i),A3(t)) , (10) 

with C/t (e e(^\ e^'^)), where e*^*' denotes the z-th 
eigenvector corresponding to the eigenvalue \i{t). In- 
serting the expression Eq. (|10|l into Eq. Q one finds 
(repeated indices are summed) 

{{uuxiit)) {u.uxm)^) ^ {ut)p^x,m 

= 5,,\dt) . (11) 

Thus the optimal operators Xi which have maximal over- 
lap with the physical states are obtained as linear com- 
binations of our original operators Xi, 

X. = Y.^^^ ^it^ 4'^=^^*- (12) 

j 

Thus the entries of the «-th eigenvector e provide the 
mixing coefficients for the optimal operator Xi coupling 
to the i-th state. 

The coefficients c'j^ are not only interesting from a 
technical point of view for constructing the optimal oper- 
ator, but also provide interesting insight into the physics 
of the nucleon system. These physical aspects will be 
addressed in detail in the next section. 

With the normalization chosen for our operators Xi the 
coefficients Cj are real within error bars, i.e. the imagi- 
nary parts are at least two orders of magnitude smaller 
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FIG. 4: The mixing coefficients Ci from the first and second 
eigenvector as a function of t. The data are for arUq = 0.12, 
16^ X 32. 



than the real parts and in the foUowing we only use the 
real parts 34] . In Fig. ^ we show the coefficients Cj de- 
termined from the first and second eigenvector which de- 
termine the optimal operator for the ground state and 
the excited state. 

Strictly speaking the matrix U which is used in Eq. 
(|10|l to diagonalize the correlation matrix M{t) is time 
dependent. However, for a proper physical state one ex- 
pects that the coefficients are independent of t within 
the time interval that corresponds to the given physical 
state. The time independence of the coefficients is well 
pronounced for the nucleon = 3 ... 17 in the top plot) 
and also for the two negative parity states (t = 22 . . . 30). 
The regions of constant Cj agree nicely with the regions 
where we find plateaus in the effective masses. Only in 
the crossing region near t = 19 and close to i = 1 for the 
positive parity states, respectively t = 32 for the nega- 
tive parity states, where higher excited states still con- 
tribute we find strong deviation from a constant behav- 
ior. Only the excited positive parity state which is seen 
in the second eigenvector effective mass plot of Fig^does 



not show long regions with t-independent coefficients. It 
seems that this channel is for small t dominated by an ex- 
cited state between t = 2 and t = 7 and then changes to 
a lower-lying state with different values for the Ci visible 
until t ~ 17. This additional mixing is clearly reflected 
in the problems which we encountered when trying to fit 
the mass of the state N'^ . The corresponding low- lying 
artifact was discussed in Section 4 and the corresponding 
structure in the correlator was shown in Fig. |2| 

An interesting question is how the coefficients depend 
on the quark mass or pion mass. In order to address 
this question we take for each state at each value of the 
quark mass the values of the corresponding coefficients 
in the center of the effective mass plateau. In particular 
we chose t = 8 for the coefficients of the nucleon, t — 28 
for the two negative parity states (this value changes to 
t = 20 on the 12^^ x 24 lattice) and t = 5 for the state 
N"^. In Fig. Owe plot the coefficients as a function of the 
pseudoscalar mass (the data for Mps were taken from 
[15|) and compare the results from 16'^ x 32 (symbols) 
to 12"^ X 24 (thin crosses) in order to check for finite size 
effects. 

The absence of finite size effects which was already 
seen for the masses is now also supported by our findings 
for the coefficients. The results from 12'^ x 24 (crosses) 
agree with the data from 16'^ x 32 within error bars. It 
is obvious that the coefficients show a non-trivial mass 
dependence. The coefficients necessary for the construc- 
tion of the optimal operators Xi which couple to the z-th 
state can be read of directly from the figure. We remark 
that the optimal operators Xi are a powerful tool which 
can be used to study properties of excited states other 
than just their mass such as e.g. matrix elements. 



VII. PHYSICAL INTERPRETATION OF THE 
MIXING COEFFICIENTS 

As mentioned, the time component of operator X3 cou- 
ples to the spin 1/2 state. In order to check its influence 
we also study the (xi, X2)-subsector alone. Fig. gives 
the resulting mixing coefficients analogous to Fig. |S1 We 
find, that xi and X2 behave similar in both cases and 
that X3 plays a subordinate role. Due to the different 
definition (e.g. our states are orthogonal and normal- 
ized) it is difficult to compare the numbers of Fig. Elwith 
the results of In both cases, as in 5], the nucleon is 
overwhelmingly due to xi while N'^ is mostly due to X2- 

First, we would like to give a simple explanation, based 
on a quark model picture, why the interpolator X2 cou- 
ples only weakly to the nucleon. According to the quark 
model the nucleon belongs to the 56 representation of 
SU{6). Its flavor-spin wave function is completely sym- 
metric with respect to permutation of two constituent 
quarks. Since the color wave function is completely an- 
tisymmetric, the Pauli principle requires that the wave 
function must by symmetric with respect to permutation 
of the spatial coordinates of any two constituent quarks. 
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Hence the parity of any diquark subsystem is positive. 
Consequently the pure SU{6) symmetric nucleon wave 
function does not contain any pseudoscalar or vector di- 
quark subsystems and the nucleon cannot couple to the 
interpolator X2- This argument is exact in the large 
limit. 

The other question is why the nucleon couples to xi 
much stronger than to X3 i^i spite of both having simi- 
lar diquark content. The answer is that the interpolators 
belong to different chiral representations (as discussed 
in Sect. 3). If chiral symmetry was not broken then each 
hadronic state would couple to the interpolator with con- 
crete transformation properties under SU{2)l x SU{2)b,- 
Chiral symmetry breaking implies that the nucleon does 
not belong to any specific representation and couples 
to X1jX3 and weakly also to X2, all of them belong- 
ing to distinct representations. The stronger coupling 
to xi indicates that the corresponding representation 
(0, 1/2) -t- (1/2,0) dominates the nucleon wave function. 

The SU (6) symmetry is broken by the residual spin- 
spin and other interactions, which are corrections of 
0{1/Nc). The weight of the components that contain 
pseudoscalar or vector diquark subsystems is very small 
(e.g. for a chiral quark model cf. Ref. [2^). This explains 
why the component C2 is small in the panel for the nu- 
cleon in Fig.[Sland why one does not see the nucleon from 
the interpolator X2- While we do observe some coupling 



of the nucleon to the interpolator ^3 the corresponding 
amplitude C3 in the upper Ihs. plot is small. We stress 
that the Ci are amplitudes, i.e. their squares are proba- 
bilities adding up to 1. 

Within the quark model the Roper state belongs to 
the other 56 representation of SU{6). The Roper is a 
radial excitation of the spatial part of the nucleon wave 
function. Hence all the arguments given above for the 
nucleon apply also to the Roper. However, the coeffi- 
cients plot for the N'_^ (lower Ihs. plot in Fig. [Sj) shows a 
large component C2 and assuming the quark model pic- 
ture one concludes that the N'^ state seen with either 
the diagonal X2-correlator or with the eigenvalue A2 is 
not the Roper state. A possible candidate for the N'^ 
state could be some higher-lying positive parity state as 
e.g. the A^(1710), which does contain a large pseudoscalar 
diquark component in its wave function. 

Another important issue which we would like to discuss 
is the pronounced quark mass dependence of the mixing 
coefficients Ci for the iV(1535) and A^(1650) states, shown 
in the two rhs. plots of Fig. [S] One clearly sees that 
while at large quark masses the A^(1535) is dominated 
by xi and the A^(1650) by X2, the situation becomes 
different towards the chiral limit. This feature persists 
also when studied just a subset of operators as in Fig. 
El Upon extrapolation to the chiral limit the A^(1535) 
couples optimally (up to some uncertainty) to (xi ~ X2), 
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while 7V(1650) couples to (xi +X2) and the contribution 
of X3 is suppressed for both negative parity states. The 
very fact that the mixing towards the chiral limit is quite 
different from that in the heavy quark region implies that 
the physics of these resonances (and in particular of their 
wave functions) is strongly influenced by chiral symmetry 
and its spontaneous breaking. It is interesting to try to 
connect this behavior to the quark model. Both A''(1535) 
and iV(1650) belong to the negative parity L = 1 70-plet 
of SU{6). Each member of this multiplet contains in its 
wave function both scalar and pseudoscalar diquark com- 
ponents. The mixing of these two different components 
is provided by tensor forces. 

A proper mixing of two independent components does 
explain, in particular, peculiarities of the N* — > Nrj de- 
cays. One can construct a linear combination such that 
the Nr] decay of the A^(1650) gets strongly suppressed, 
while the decay of the N{1535) gets enhanced. There 
are two possible origins of a tensor force providing the 
necessary mixing: The tensor force which is supplied by 
perturbative gluon exchange and which does not contain 
flavor dependence in its operator structure [2^ and the 
tensor force mediated by the Goldstone field and related 
to spontaneous chiral symmetry breaking 0, |23| [l^ . If 
the mixing was due to perturbative gluon exchange one 
would not expect that the mixing would depend on the 
current quark mass. However, a tensor force that orig- 
inates from chiral symmetry breaking appears only to- 



wards the chiral limit. Our lattice results, shown in Figs, 
m and El clearly support the latter possibility, though 
these results cannot shed any light on the specific micro- 
scopical origin of the tensor force. This result is also con- 
sistent with the recent phcnomenological '2^ and large 
Nc 29] analyses of the mixing of the negative parity 
baryons: Both studies suggest that a proper mixing is 
supplied by the flavor-dependent tensor force. 



Finally, our results allow to explain the smallness of 
the 7riViV(1535) and 7riViV(1650) coupling constants as 
compared to the large ttNN coupling. It has been shown 
in Ref. [s^ that if both the nucleon and A^(1535) (or 
iV(1650)) are created from the vacuum only by the xi 
and X2 interpolators, irrespective of their mixing, then 
the chiral symmetry properties of these interpolators im- 
ply that the 7r7V7V(1535) and 7riViV(1650) coupHngs must 
vanish. In the present work we find that the coupling 
of all these baryons to the interpolator X3: which be- 
longs to another representation of the chiral group, is 
much smaller than the coupling to xi or X2- This qual- 
itatively explains the smallness of the 7rA^A^(1535) and 
7rA^A^(1650) coupling constants. 
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VIII. CONCLUSIONS 

Let us summarize our main findings. We use three dif- 
ferent interpolating fields and study how these fields cou- 
ple to the nucleon and its excited states of both parities. 
We clearly identify distinct negative parity excited states 
iV(1535) and A^(1650) and show that the mixing of dif- 
ferent interpolators which create these baryons from the 
vacuum, is obviously quark mass dependent. The mixing 
towards the chiral limit is qualitatively different from the 
mixing in the heavy quark region. This implies that the 
physics of these states is strongly related to chiral sym- 
metry of QCD. A very specific mixing which is observed 
towards the chiral limit can be a basis for the explanation 
of N* Nt] and N* Ntt decays of these baryons. 

In a quark model picture we are able to understand 
peculiarities of the coupling of different interpolators to 
the nucleon and its excited states. 

We clearly see an excited state of positive parity at 
relatively large quark masses from the interpolator X2- 
Although we cannot reliably trace this state towards the 
chiral limit we argue that it lies too high for an identi- 
fication as the Roper resonance and may correspond to 
the A^(1710) state. 

While we do observe some contribution of positive par- 
ity excited state(s) to the xi interpolator at small times, 
the signal is rather weak and using traditional meth- 
ods, such as effective mass plots or one (two, many)- 
exponential fits of the correlator, do not allow us to con- 
clude whether this signal comes from the Roper state or 
from some higher lying states. It would be an important 



task to optimize source and sink in order to improve the 
signal from the Roper resonance. 

In the subtraction of a parametrization of possible 
quenched artifacts [r]' N) — giving rise to negative val- 
ues of the correlation function for central ^-values — was 
crucial for the identification of a Roper signal. We also 
observe such negative values but postpone an analysis to 
future studies with better statistics. 

An interesting challenge for further lattice studies is 
the empirically observed approximate parity doubling of 
nucleons above 1.7 GeV region (and similar in delta spec- 
trum). It has been suggested that this doubling may 
reflect a smooth effective chiral symmetry restoration in 
the upper part of hadron spectra 23] . Still, there remains 
the possibility that this doubling is accidental. Lattice 
methods can help to clarify this interesting issue. 
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